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ABSTRACT 

The nearby, middle-aged PSRB1055—52 has many properties in common with the Geminga pulsar. 
Motivated by the Geminga’s enigmatic and prominent pulsar wind nebula (PWN), we searched for 
extended emission around PSRB1055—52 with Chandra ACIS. For an energy range 0.3 — IkeV, we 
found a 4 a flux enhancement in a 4"9 — 20" annulus around the pulsar. There is a slight asym¬ 
metry in the emission close, 1"5 — 4", to the pulsar. The excess emission has a luminosity of about 
10 29 ergs -1 in an energy range 0.3 — 8 keV for a distance of 350pc. Overall, the faint extended emission 
around PSRB1055—52 is consistent with a PWN of an aligned rotator moving away from us along 
the line of sight with supersonic velocity, but a contribution from a dust scattering halo cannot be 
excluded.Comparing the properties of other nearby, middle-aged pulsars, we suggest that the geome¬ 
try - the orientations of rotation axis, magnetic field axis, and the sight-line - is the deciding factor 
for a pulsar to show a prominent PWN. We also report on an > 30% flux decrease of PSRB1055—52 
between the 2000 XMM-1V ewton and our 2012 Chandra observation. We tentatively attribute this 
flux decrease to a cross-calibration problem, but further investigations of the pulsar are required to 
exclude actual intrinsic flux changes. 

Subject headings: pulsars: individual (PSRB1055—52, Geminga) — stars: neutron 


1 . INTRODUCTION 

When pulsars spin down, most of their rotational 
energy loss is carried away by winds of relativistic 
particles. When the winds shock in the ambient 
medium, these shocks emit synchrotron radiation which 
becomes observable as a p ulsar win d nebula ( PWN, 
for re views see, e.g., I Slang 1201 It IKarealtsev &: Pavlov! 
120081b Gaensler fc Sland I2006l i7 The PWN luminosity 
Lpwn correlates with the pulsar spin-down power 
E as Lpwn = E x 7ypwN- The efficiencies t?pwn at 
X-ray ener gies are usually betwe en rjx ~ 10 -5 and 
? 7 x ~ 10 ~ 3 iKargaltsev et all [2012 1. The wide range of 
observed efficiencies depend on ages, magnetic fields, 
environments, and evolutionary states of the neutron 
stars. Older pulsars have smaller E, their PWNe are 
dimmer and more difficult to observe. Thus, most of 
the currently known ~ 60 PWNe are powered by young 
(< lOOkyrs) pulsars, while little is known about PWNe 
around older pulsars. Only very few PWNe around 
middle-aged (0.1 — lMyr) pulsars with E < 10 35 ergs -1 
have been detected. 

The best known nearby representatives of such 
middle-aged pulsars are PSRB1055—52, Geminga and 
PSR B0656+14. They have similar spin-down energies 
(E ~ 10 34 ergs -1 ), similar ages (a few lOOkyrs), and 
similar magnetic fields (~ 10 12 G). All three are 7 -ray 
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pulsars. Their X-ray spectra consist of non-thermal 
and strong thermal components. Because of the sim- 
ilarities, they were dubb ed the Three Musketeers by 
iBecker fc Triimoerl ( 1997 1. The X-ray spectra are 
best described by a three-component model - a cold 


and a hot blackbody plus a power law componen t 
(e.g.. IBecker &. Aschenbachl l2002t iPavlov et all ) 20021) . 
iDe Luca et all ( 2005f) studied the individual components 
with phase-resolved, low-resolution X-ray spectroscopy. 
They found pronounced differences for these three 
sources indicating different orientations, magnetic field 
configurations and surface temperature nonuniformities. 
The current (June 2015) statistic of nearby (< 2kpc), 
middle-aged pulsars (characteristic ages between 
0.1 — lMyr and E between 5 x 10 33 — 5 x 10 34 ergs -1 ) 
lists 17 objects (see Table [U, many recently discovered 
with Fermi. 


Putative PWNe around the Three Mus keteers , 
found with ASCA (e.g., iShibata et all 1997 1. were 
later resolved into several point sources, thanks to the 
b etter spatia l reso lution of ROSAT and BeppoSAX 
(jBecker et al.l I1999H . Further improvement in spatial 
resolution and sensitivity resulted in the XMM-Newton 
detection of two long , lateral PWN-tails behind Geminga 
(jGaraveo et all 120031 1. The superb spatial resolution of 
Chandra revealed the axial tail and the structure of the 
lateral tails in the immediate v icinity of the Geminga 
pulsar (|Pavlov et al.ll2010ll2006ll . Of the 12 middle-aged 
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pulsars (except for PSRB1055—52) with sensitive, 
public X-ray observations, 6 have prominent X-ray 
PWNe, while the other 6 have no detected or very faint 
PWNe; see Table Q] PSR J1741-2054, for example, 
is surrounded by a large H Q bo w shoc k ne bula an d 
an impressive X-ray PW N dBrownsberger fc Romanil 
l2?TIl IRomani et al.ll2?)Ioh . PSR J0357+3205 has a very 
long tail, though it was ruled out that this long tail 
is associated with a classical bo w-sh ock, ram-pressure- 
dominated PWN (e.g., iMarelli et alJ 120131 ). In the 
case of PSR B0656+14, Chandra observations found 
evidence for a very faint or compact PWN (Birzan et 
al. 2015, submitted). No significant extended emis¬ 
sion was seen for PSRB1055—52 with XMM-Newton 
MOS 1/2 dDe Luca et al.l 120051 : iBecker fc Aschenbachl 
2002). The Contin uous Clockin g mode observation of 
PSRB1055—52 by iPavlov et all (|2002t ) did not allow 
for exploiting the full spatial resolution capabilities of 
Chandra to search for faint extended emission. Our new 
Chandra ACIS observation, reported here, enables such 
an investigation. 

We want to note that the distance estimate for 
PSRB1055—52 has been revised and is dif feren t from the 
value reported in the ATNF catalog (Mancheste r et all 
I2005H . The dispersion measure of PSRB1055—52 con¬ 
verts to either d = 1.53 kpc using the Galactic free 
electron density model by iTavlor fo Cordesl ( 1993 ). or 
to d = 0.73 ± 0.15k p c using the NE2001 model by 
iCordes fc Laziol (12002D . iMignani et, al.l (|201flH analyzed 
optical emission from the pulsar together with its X- 
ray spectrum. Considering the contributions of the 
individual thermal components to the multiwavelength 
spectrum and an upper limit on the neutron star ra¬ 
dius of 20 km, they inferred a much smaller distance, 
350 ± 150 pc, than the DM-based ones are. Thus, 
PSRB1055—52 is very likely comparably close to us as 
the other Musketeers are. Above all, this distance value 
also solves the previous puzzle of an unusually high ef¬ 
ficiency at y-ray waveleng ths. Scaling the Fermi LAT 
fluxes reported bv lAbdo et all (12010 1 with the new dis¬ 
tance shows that the y-ray efficiency of PSRB1055—52 
is not so different from other y-ray pulsars. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Chandra 

PSR B1055—52 was observed with Chandra ACIS-I in 
imaging mode in VFAINT telemetry mode on 2012 Octo¬ 
ber 3rd for 56 ks (ObsID 13789). We used CIAO (version 
4.6.1.) with CALDB (version 4.5.9) for the data reduc¬ 
tion and analysis. 

We investigated the possibility of pile-up0 since this 
could influence the spectral analysis and the input 
for the point spread function (PSF) modeling with 
MARX0. Using the CIAO tool pileup_map, we obtained 
0.09 photons per frame of 3.2 s in the center of the pulsar 
image (binned to 1 native ACIS pixel, 0'/49). This 
number translates into a rough estimate of the amount 
of pileup fraction, 4%, according to the approach by 

1 www.atnf.csiro.au / research / pulsar / psrcat 

2 Two or more photons are detected as a single event; for more 
details see cxc.harvard.edu/ciao/ahelp/acis_pileup.html 

3 space.mit.edu/cxc/MARX/index.html 


iDavisl (120011) . This amount of pileup is unlikely to 
significantly influence the pulsar spectrum. Hence, we 
neglect pile-up for our further data analysis. 

The data were free of background flares. Using 
specextract, we obtained the pulsar spectrum from a 
circular region with radius 5 ". We binned the data with 
each bin containing at least 25 counts. 

In order to evaluate our data for extended emis¬ 
sion aroun d the p ulsar, we used the ray-trace simulator 
MARX (iDavis et al.l [2012 ) to model the Chandra PSF. 
MARX (version 5.0) and the asso ciated Chandra Ray 
Tracer (ChaRT: iCarter et al.l 120031) program require an 
input spectrum for accurate modeling of the pulsar PSF. 
The MARX calibration data are based on CALDB 4.4.7. 
Since the contamination on the optical-blocking filters 
of the ACIS detector^! is changing the response of the 
instrument over time, mixing of the CALDB versions in 
the data and MARX simulations would produce inaccu¬ 
rate results. Therefore we extracted the pulsar spectrum 
for the MARX modeling from data re-processed with 
CALDB 4.4.7. The spectral model files were created 
from spectral fits using Sherpa and XSPEC. Following 
standard CIAO threads, we tested different approaches. 
In one, the output of ChaRT is used in combination with 
MARX. In another, the simulation is done solely within 
MARX. We obtained significantly more counts in the 
central region of the pulsar using ChaRT+MARX than 
we measured in the actual data. In contrast, the MARX- 
only simulation produced count numbers comparable to 
the observed ones. Therefore, we only discuss the latter 
simulations in this paper. In order to minimize the sta¬ 
tistical errors of the simulated PSF, we carried out 1000 
MARX simulations with an exposure time of 1 Ms using 
the internal MARX dither model. Our final MARX pul¬ 
sar PSF model is the average of these 1000 simulations, 
scaled to the actual exposure time (56ks), the consid¬ 
ered simulation ‘errors’ are the count errors scaled to the 
actual exposure time as well. 

2.2. XMM-Newton 

Our spectral fit results for the Chandra pulsar data 
deviate from previously published results (Section 13.21) . 
In order to directly compare our Chandra and previous 
XMM-Newton data, we reanalyzed the XMM-Newton 
observations (ObsIDs 0113 050101, 0113050201) from 
2000 which were presented bv De Luca et al.l (|2005h . The 
EPIC (European Photon Imaging Camera) observations 
of PSRB1055—52 employed the pn camera in timing 
mode, and MOS1 and MOS2 cameras in imaging mode. 
For all instruments the medium filter was used. We used 
the XMM-Newton Science Analysis Software (SAS), ver¬ 
sion 13.5 for the re-processing of the observations and 
data reduction. These early XMM-Newton EPIC-pn 
observations do not have associated offset maps, and 
offset map correction or soft photon noise reduction 
is not possible. Soft background n oise can influence 
the extracted spectr a below 400 eV dBurwitz et al.l[20041 : 
iDe Luca et al.ll2005l) . hence the extracted pn spectra only 
cover energies larger than 400 eV. 

4 See e.g., cxc.harvard.edu/ciao/why/acisqecontam.html 
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TABLE 1 

Nearby middle-aged pulsars 


Name 

L/Q a 

Age 

kyr 

E33 h 

Bn b 

Dist 

kpc 

Ref 3 

PM 

masyr -1 

ReF 

^trans 

kms -1 

N h 

10 20 cm - 2 

ReF 

r d 

1 7 

E m t d 

GeV 

PWN e 

J0357+3205 

RQ 

540 

6 

24 

0.6 

(gd) 

164 

( 6 ) 

466 

8 ± 4 

(13) 

1.0 ± 0.1 

0.8 ± 0.1 

Y 

J0358+5413 

RL 

564 

45 

8 

1 

(P,2) 

12.3 

(7) 

58 

< 72 

LAB 



Y 

J0538+2817 

RL 

618 

49 

7 

1.3 

(P,2) 

57.89 

( 8 ) 

357 

< 54 

LAB 



Y 

J0622+3749 

RQ 

208 

27 

30 

1.9 

(gd) 




< 21 

LAB 

0.6 ±0.4 

0.6 ± 0.1 

N 

J0633+1746 

RQ 

342 

32 

16 

0.25 

(P,2) 

169 

(9) 

200 

1.07 

(14) 

1.2 ± 0.1 

2.2 ± 0.1 

Y 

J0659+1414 

RL 

111 

38 

47 

0.28 

(P>2) 

44.1 

( 10 ) 

59 

4.3 ±0.2 

(14) 

1.7 ±0.5 

0.4 ±0.2 

VF 

J0922+0638 

RL 

497 

7 

25 

1.1 

(p> 2 ) 

88.4 

( 11 ) 

461 

< 3.2 

LAB 



N 

J1057—5226 

RL 

535 

30 

11 

0.35 

(3) 

42 

(3) 

70 

2.7 ±0.2 

(14) 

1.1 ± 0.1 

1.4 ±0.1 

VF 

J1741—2054 

RL 

386 

10 

27 

0.38 

(P,4) 

109 

(4) 

196 

12 . 0 +°? 

(4) 

1.1 ± 0.1 

0.9 ±0.1 

Y 

J1745—3040 

RL 

546 

9 

20 

0.2 

(p> 2 ) 

7 

( 12 ) 

7 

< 99 

LAB 



N 

J1846+0919 

RQ 

360 

34 

15 

1.5 

(g,l) 




< 33 

LAB 

0.7 ±0.3 

2.2 ±0.5 

N 

J2028+3332 

RQ 

576 

35 

9 

0.9 

(g.5) 




< 62 

LAB 

1.2 ± 0.2 

1.9 ±0.3 

N 

J2030+4415 

RQ 

555 

22 

12 

0.8 

(g.5) 




6 + 15 

(5) 

1.6 ± 0.1 

1.7 ±0.3 

Y 

J1530—5327 

RL 

944 

9 

12 

1.46 

DM 




< 65 

LAB 



? 

J1549—4848 

RL 

324 

23 

20 

1.54 

DM 




< 44 

LAB 



?(ND) 

J1746—3239 

RQ 

482 

33 

12 






< 55 

LAB 

1.4 ±0.1 

1.5 ±0.2 

? 

J1957+5033 

RQ 

839 

5 

17 

0.9 

(gd) 




< 12 

LAB 

1.3 ±0.2 

1.0 ± 0.2 

?(ND) 


a Radio-quiet (RQ) or radio-loud (RL); based on the ATNF catalog fwww.atnf.csiro.au/research/pulsar/psrcat; IManchester et alj|2005 
b The spin-down energy, E 33 , is in units of 10 33 ergs -1 . The magnetic field Bn is in units 10 11 G. 

c References: (g,l) indicates that the distance was estimated from 7 -ray LAT flux as described by (1) Brownsberger &; Romani (2014); 
(p,2) denotes parallactic distances corrected for the L utz-Kelker bias as listed by (2) IVerbiest et a l. (2012); (3) Mignan^et al. ffOlO); DM 
dis tance estimate based on the free electro n mod el by 
(5) IMar elli et ahl (120150 ; ( 6 ) IDe Luca et I 
IBrisken et al.l (I2003bl) ; (11) [B r 

the reference ‘LAB’ indicate upper limits which are the total Galactic HI colu mn density in the direction of the pulsar using the nh-webtool 
(heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl) and the LAB survey dKalberla et al.1120051) 

d The spectrum of each Fermi LAT pulsar was modeled as a power law with an exponential cutoff by lAbdo et 31 mm . their ‘PLECr 



model. In this table we quote their reported photon index, r ELEC1 , and the cutoff energy, E^ EC1 . The cutoff was significantly detected for 


all pulsars - all have reported cutoff significance > 33. 

e Y - prominent PWN; N - no PWN detected in current XM.M-Newton or Chandra data; VF - indication of very faint PWN; ? - unknown 
(non-existent X-ray observations with XMM -Newton or Chandra ); (ND) - new data expected from scheduled/recent X-ray observations 


20 " 



V.'' E- 

Fig. 1.— Our analysis regions around PSRB1055—52 which we 
used to check for large-scale extended emission. The circle on the 
top marks the maximum radius ( 20 // ) we used for annuli around the 
pulsar. For these measurements background region 1 (BG1) was 
applied. The fan structure targets larger radii from the pulsar, but 
avoids contaminating sources A and B, as well as ACIS chip gaps. 
The image is smoothed and covers the energy range 0.3 — 8 keV. 


We checked the Chandra image for nearby X-ray 


sources that could contami nate t h e XMM-Newto n data 
(see Figure |I}. As noted bv IDe Luca et al.l (! 2005l ). there 
is a faint source (source A) in northeast direction. Its 
angular separation from the pulsar position is 30" in 
our Chandra observation. In the pn timing mode it is 
< 4 pixels along the RAWX direction, whe re spati al 
inform ation is maintained for this mode (IDe Luca et al.1 
120051 1. In Section 13.21 we use the Chandra data to assess 
the influence of this contaminating X-ray source on the 
( XMM-Newton ) pulsar spectral parameters. There is 
also a fainter source B at 39" east of the pulsar which 
is not distinguishable from the broad pulsar PSF in the 
MOS image. There is furthermore a very faint source 
C in the northeast direction with an angular separation 
from the pulsar of 3', but it is projected on a similar 
RAWX pixel as source A in the pn timing observation. 
Source C is a factor 4 fainter than source A, though. For 
the spectral extraction in the pn timing observations 
we used 33 <RAWX< 39 (both observations 101, 201) 
for the source regions, and 5 <RAWX< 7 (101) and 
4 < RAWX < 6 (201) for the background regions. For 
the MOS cameras, the source region is a circle centered 
on the pulsar with radius 45", excluding two circular 
regions (radius 10") centered on sources A and B. 
The background spectra were extracted from nearby 
source-free regions on the same detector chip. After 
filtering for background flares, we extracted spectra 
using events with pattern 0 to 12 for MOS, but only 
single and double photon events (pattern 0 to 4) for 
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Fig. 2.— PSRB1055—52 and its close surroundings in X-rays 
(0.3 — 8keV) as seen by Chandra. The arrow points in the 
direc tion of proper motion 94° ±7° east of north lIMignani et al.l 
l201 0f). The annulus region has a size of 1.5” — 4 ,/ . The four 
quadrants used for count measurements are labeled. The small 
wedge inside the annulus indicates the Chandra PSF asymmetry 
region within 1". 


EPIC-pn. The XMM-Newton spectra were binned using 
at least 40 counts per bin. 


2.3. Spectral analysis 

Spectral fits were done using XSPEC version 12.8.1. 
We excluded energies higher than 7 keV since the 
background dominates at these energies. We used 
the Tuebingen-Boulder ISM absorption m odel (tbabs) 
with the solar abundance table from 1 Wilms et all 
(l2000f). the photoelectric cross - sectio n table from 
iBalucinska-Church fc McCammonl (119921) toget her with 
the He cross-section based on lYan et al.1 (119981) . Note 
that instead of using merged XMM-Newton data sets, we 
obtained simultaneous spectral fits for all the individual 
six XMM-Newton spectra. 


3. RESULTS 

In the Chandra ACIS data, in a circular region cen¬ 
tered on the pulsar with radius 2" and 5" we detected 
1858 ±43 and 1908 ±44 (background subtracted) counts, 
respectively, in the energy range 0.3 — 8 keV. There is no 
obvious large scale extended emission (Figure [TJ but a 
slight asymmetry of the count distribution within 4" of 
the pulsar (Figure [2]). 

3.1. Extended emission around the pulsar 

We consider three different energy bands for our ex¬ 
tended emission analysis: the broad energy band (0.3 — 
8keV, BEB), the soft energy band (0.3—1 keV, SEB), and 
the hard energy band (1 — 8keV, HEB). As explained in 
Section [21 we employ for comparison MARX simulations 
for 1 Ms, scaled to exposure time of the observation. We 
tested different ‘AspectBlur’ parameter values and found 
0'/2 to best reproduce the radial count distribution of the 
inner PSF of the observed pulsar. 


TABLE 2 

Significance of the radial flux enhancement. 


0.3 - 

8.0 keV 

0.3 - 

1.0 keV 

1 - 8.0 keV 

Obs 

Sim 

Obs 

Sim 

Obs Sim 


rin [arcsec] 

4.88 

4.61 

4.88 

4.61 

4.88 

4.61 

r out [arcsec] 

13.19 

14.21 

20.00 

20.00 

10.01 

10.01 

N ( a ) 

73.00 

10.30 

39.00 

6.08 

30.00 

3.55 

SN 

8.54 

3.21 

6.24 

2.47 

5.48 

1.88 

A [arcsec 2 ] 

472 

568 

1182 

1190 

240 

248 

n 0>> 

0.16 

0.089 

0.033 

0.010 

0.13 

0.081 

5n 

0.02 

0.008 

0.005 

0.002 

0.02 

0.008 

^obs ^sim 

0.066 

0.0231 

0.044 

^(^obs ^sim) 

0.020 

0.0057 

0.024 

S< c > [a] 

3.3 

4.1 

1.8 


Source counts for the observation (“Obs”) and scaled counts 
of the 1 Ms MARX simulation (“Sim”) 

C) Surface brightness [cts arcsec -2 ] includes background counts 
for the observational data. The background surface brightness 
is added to the scaled surface brightness of the 1 Ms MARX 
simulation 

Surface brightnesses, n, their differences, uncertainties, and 
the significance, S, were calculated before rounding the respec¬ 
tive contributing values. 

3.1.1. Radial enhancement of surface brightness 

The radial distributions of the surface brightness 
of the MARX simulation and the observation are 
shown in Figure [3] For this plot, the background is 
not subtracted from the observational data. Instead, 
the respective background surface brightness of the 
observation was added to the MARX simulation. The 
surface brightness of the background is 0.071, 0.005, 
and 0.066 counts arcsec -2 for the BEB, SEB and HEB, 
respectively. The uncertainty leveQ of the MARX 
simulation curve at the largest distances is actually 
dominated by the uncertainty of the added observational 
background surface brightness. 

There is an enhancement of X-ray emission starting 
from about 2" to about 20". It is especially well visible 
in the soft band. We evaluated the significance of the 
enhancement in Table [2] Since the MARX simulation 
allowed a finer binning, the considered observation and 
simulation annuli do not have exactly the same size. For 
the comparison in Table [2j we used for the simulation 
annuli with areas at least as large or larger than the 
areas of the annuli utilized for the observational data. 
This ensures that our significance estimate in Table [2] 
represents a lower limit to the actual significance of the 
enhancement in the considered region. The enhancement 
in the surface brightness has a significance of 3 a in the 
BEB and 4<r in the SEB, while in the HEB the surface 
brightness in the observation is consistent with the 
simulation. We checked that the enhancement persisted 
even if we used different (and badly fitting) ‘AspectBlur’ 
values (the probed range was 0'/07 — 0"8). The reason 
is that only the inner part (within « 4") of the PSF is 
significantly changed by this MARX parameter. 

5 The error for the simulated surface brightness, n s i m , was calcu¬ 
lated as: <5n sim = A~ 1 (N^P n t exp /t siln + AT b g 41 2 /A| g ) 1 / 2 j where 
texp is the exposure time, t s [ m is the simulation exposure time, 
is the number of counts scaled by the exposure time, TVbg i s 
the number of background counts, A is the area of the respective 
annulus, and Abg is the area of the background region. 
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Fig. 3. — The radial surface brightness distribution in the obser¬ 
vation (blue) and the MARX simulation (red; 1 Ms exposure time, 
scaled to source exposure time, internal dither model, AspectBlur= 
O'/2) in the broad energy band (top panel), the soft energy band 
(middle panel), and the hard energy band (bottom panel). The 
‘error’ bars in distance direction indicates the bin size used for the 
surface brightness measurements, the error bar in y direction is 
the actual uncertainty for the respective bin. Note that no back¬ 
ground subtraction was done for the observational data. Instead, 
the respective background level was added to the MARX simula¬ 
tion surface brightness. 


In order to crudely estimate the flux in the X-ray 
emission enhancement in the BEB, we applied the 
ef f 2evt tool which uses the event energies and effective 
areas at the event locations for the flux estimate. We 
used the 4.88" — 13.19" annulus (as in Tabled]) for the 
data and the MARX simulation, and did the same for 
the (area-scaled) background region. The resulting flux 
for the excess emission is 1.0 ± 0.2 x 10 -14 ergs -1 cm' 2 
for the BEB (0.3 — 8keV). Alternatively, using the 
count rate of the BEB enhancement in the same an¬ 


nulus (5.5 x 10 4 cps), PIMMS0 with Chandra’s A013 
response, and an absorbed (Ah = 3.4 x 10 2O cm -2 
from the XMM /Chandra fit of the pulsar spectrum in 
Table 0]) power law with photon index range of T = 2 
and F = 4, we obtain absorbed (unabsorbed) BEB fluxes 
of 0.6 x 10” 14 (0.7 x 10 -14 ) erg s _1 cm -2 and 1.1 x 10~ 14 
(1.7 x 10 -14 ) ergs -1 cm -2 . 

We checked whether there is any flux enhancement at 
angular separations 20" to 60". For this we use fan¬ 
shaped regions to avoid chip gaps. Although we used 
a different background region (see Figure Q]), the back¬ 
ground rates from this region are consistent with the ones 
quoted above for the BEB, SEB, and HEB. We did not 
find any significant flux enhancement for angular sepa¬ 
rations > 20". 

3.1.2. Probing for asymmetry in the extended emission 

Figure [2] appears to show a slight asymmetry of the 
count distribution within 4" of the pulsar. This asym¬ 
metry is notably in the direction of the PSRB1055—52 
proper motion. Aiming to investigate this asymmetry, 
we compare in Figure Q] the East and West side of the 
close neighborhood of the pulsar (> 1"5 to exclude the 
known PSF asymmetry regiorQ). There are more counts 
in the second annulus (center at 2"6) for the East data 
than for the West data. Count numbers for these 10 
annuli are, however, so small that 2 a errors overlap for 
the annulus in question. Thus, the difference between 
East and West is not statistically significant. 

Based on our MARX simulation, we also tried an¬ 
other approach. We divided the immediate surrounding 
of the pulsar in four annular (1"5 — 4") quadrants (Ql: 
319°—49°, where 0° and 90° indicate the North and East 
direction, respectively, Q2: 49° — 139°, Q3: 139° — 229°, 
Q4: 229° — 319°; see Figure [2]). Similarly to the an¬ 
nuli estimates in Table [5] in Section 13.1.11 we calculated 
the difference of the surface brightness considering the 
observed and the simulated MARX data in the BEB, 
SEB, and HEB in each quadrant. The observed respec¬ 
tive background surface brightness of the observation was 
added to the simulation data for the comparison. While 
we indeed find a slight enhancement in the second quad¬ 
rant as expected from Figure [5] the significance is just 
about 3er for the BEB and SEB for the available small 
number of counts (24 counts maximum in Q2 of which 10 

Radial Profile 10 Annuli Broadband 
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o 
U 
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Fig. 4.— Counts in the East and West direction over distance 
from the pulsar in the BEB. The measurement annuli start at l'/5. 

6 cxc.harvard.edu/toolkit/pimms.jsp 

7 see cxc.harvard.edu/ciao4.6/caveats/psf_artifact.html 
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TABLE 3 

Significance of extended emission in quadrants. 



Q1 

Q2 

Q3 

Q4 

BEB TTofog 

1.30±0.35 

2.22±0.45 

1.02±0.31 

1.39±0.36 

BEB n sim O 

1.00±0.07 

0.90±0.07 

0.73±0.06 

1.0±0.07 

BEB 

0.8 

2.9 

0.9 

1.1 

SEB T'lobs 

0.28±0.16 

1.30±0.35 

0.65±0.25 

1.02±0.31 

SEB n sim ( a ) 

0.44±0.05 

0.38±0.04 

0.30±0.04 

0.39± 0.05 

SEB S< a > 

1.0 

2.6 

1.4 

2.0 

HEB TT-obs 

1.01±0.31 

0.93±0.29 

0.37±0.19 

0.37±0.19 

HEB n sim ( a > 

0.56±0.05 

0.52±0.05 

0.43±0.04 

0.60±0.05 

HEB SO) 

1.4 

1.4 

0.3 

1.2 


( a ) The surface brightness of the simulation includes the observa¬ 
tional background component as described in the text. 
d 5 ) Significances, <S, of the surface brightness difference were cal¬ 
culated before rounding of the surface brightness values. 

are expected to be pulsar and background counts. Thus, 
we cannot firmly exclude that the hint of asymmetry for 
the emission around the pulsar in the 1"5 — 4" region is 
just a statistical fluctuation in the spatial count distri¬ 
bution. 

3.2. The pulsar 

Details on the PSRB1055—52 spectrum were pre¬ 
viously reported for the sensitive XMM-Newton 
observations (~ 120 ,500 counts for all EPIC detectors; 
iDe Luca et all [2005h . Our ACIS pulsar spectrum has 
much fewer counts (ss 1900). Fittin g the data with a 
similar spectral model as used bv lDe Luca et al.l (120051) . 
two blackbodies (BBs) and a power law (PL) with 
interstellar absorption, we obtained an acceptable fit to 
our ACIS spectrum (last column in Table 0 . We found 
differe nces with the values reported by IDe Luca et al.l 
(120051) . The difference in the absorbing column density 
is a factor 10, which cannot be explained by different 
abundances or our large uncertainties. We also found an 
indication for a larger (i.e., softer) PL index, although 
our u ncertainty range and the one of IDe Luca et al.1 
(I2005D still overlap for this parameter. In contrast 
to XMM-Newton , Chandra is able to clearly resolve 
neighboring sources. In particular, we find not only 
source A, but also an additional fainter source B located 
within the MOS PSF (Sources A arid B are labeled in 
Figured]). While IDe Luca et al.l (I2005D exclude regions 
encompassing the source A from the MOS data, this 
source cannot be completely excluded from the pn tim¬ 
ing data. It is interesting to check the influence source 
A would have on the pulsar spectrum if i t were included 
in the spectral extraction region. As IDe Luca et al.1 
(|2005l ) noted, Source A is much fainter than the pulsar, 
and it is a hard source. In our ACIS data, Source A has 
77 counts in comparison to the pulsar which has 1912 
(using circular extraction regions with radii of 5" for 
both). Using ‘cstat’ statistics and a simple PL, we find 
a photon index T = 1.5+0'? for Source A. If we fit the 
pulsar spectrum including the counts of Source A with 
the 2BB+PL model, we find a PL index T = 2.01 q 7 
versus T = 2.3+o'g if only the pulsar spectrum is fit. 
While this is not a significant change, the numbers seem 
to support the natural expectation for this hard source 
that the PL index would only get smaller if Source A 


cannot be properly removed in the spectral extraction 
region. We, however, found the opposite - the (Source 
A-free) ACIS spectral model fit gave a larger photon 
in dex tha_n what_ was found with XMM-Newton by 
IDe Luca et al.l (I2005D . 

Guided by the new spatial information from Chandra 
we reanalyzed the XMM-Newton data as described in 
Section 12.21 Due to our more rigorous source ex trac tion 
region s, we have fewer counts than IDe Luca et al.1 
( 2005 ). Overall, we obtained for our source extraction 
regions 30,644 (source and background) counts, 10,965 
counts, and 11,734 counts for p n, M OS1 and MOS2, 
respectively. As noted in Section 12.31 we did simultane¬ 
ous spectral fits for all the individual six XMM-Newton 
spectra from the two respective observations. In Ta¬ 
bled] we present our fit results for the 2BB+PL model 
considering only the pn data (second column) and a tied 
fit for the pn and MOS data. While there are differences 
in these two fits, they are usually not significant and 
well within c ross-calibration er rors of the XMM-Newton 
instruments (jR.ead et all 12014 ). The only exception is 
Nh where the pn-only fit gave a significantly higher 
value. As explained in Section [2.21 these early pn timing 
observations are a rare case where there are no usable 
counts below 400 eV in pn, but some counts for MOS. 
Since interstellar absorption has its largest effect on soft 
X-rays, Nh from the tied fit is more reliable. We find a 
slight shift to a softer P L photon index in comparison 
to the results bv IDe Luca et al.l (I2005D (1.9 ± 0.3 versus 
1.7 ± 0.1). Reasons for this small change are our more 
rigorous extraction regions, but probably also calibration 
updates for the MOS instruments. 

We fit the 2BB+PL model to the XMM-Newton data 
in combination with our ACIS data using tied Nh values 
(since ACIS is not very sensitive at the lowest energies) 
and tied photon indices (assuming ACIS might give a 
more realistic value because of better spatial resolution 
and lower background at higher energies). The results 
for these tied parameters, listed in column 4 and 5 in 
Table 2] are very close to the XMM-Newton -only fit 
(third column) as one would expect considering the 
different count statistics. An unfolded spectrum for 
the resulting best-fit model is shown in Figure [5] It is 
interesting to note that even with these two parameters 
tied, there are significant differences between ACIS and 
the XMM-Newton fit parameters. This is true for the 
BB emission areas, the cold BB temperature, but also for 
the observed flux. The observed ACIS flux in the energy 
range 0.3 — 8 keV is only 71% of the XMM-Newton EPIC 
flux in the same energy range. In Figure |H1 we show 
contour levels for the BB temperatures. Interestingly, 
even the 99% temperature contours do not overlap. We 
discuss possible reasons in Section 2^1 

As IDe Luca et al.l (120051 ) have already showed, the 
utilization of the BB model as description for the 
thermal emission of PSRB1055—52 requires at least 
two regions with different temperatures and areas and 
an additional PL component. In order to check if the 
spectrum of PSRB1055—52 can be fitted with a two- 
component (thermal+PL, or thermal+thermal) model, 
we employed NS atmosphere models for the thermal 
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TABLE 4 

Spectral fit parameters of the pulsar for different instruments. 3. 


Parameter 13 

pn Only 

pn&MOS 

pn&MOS c 

ACIS C 

ACIS Only 

N h (10 19 cm" 2 ) 

88 ± 5 

33 ± 1 

34 ± 1 

34 (tied) 

240+ 200 
Z4tu —150 

kT CBB (eV) 

61 ± 5 

Got 2 

68 l 2 

81 ± 5 

63+ 12 
° -10 

Rcbb (km) 

1 q+ 3.8 
iU -3.0 

5.3 ± 0.6 

5 .6tr 7 

9 O+0.5 

14+ 53 

±4t -13 

fclHBB (eV) 

140 ± 20 

160± 2 ° 

160j) 2 ° 

220l® 9 

170”*” 50 

11 -40 

^HBB (m) 

320” 1 ” 140 

ozu -110 

190t" 

190115 

67+35 

° * —27 

i so -1-23 ® 
J-oU-ioo 

r 

1.9 ± 0.4 

1.9 ± 0.3 

1.9 (tied) 

1.9 ± 0.2 

9 q + 0.5 

A/pl 

2 4+ 1 - 1 
z, ^-0.8 

2 Q+0.6 

Z ' U —0.5 

9 i +0.6 

Z * X -0.5 

7+0.6 

1 - '-0.7 

2 9"*” 2 ' 3 

/abs,0.3-8keV EPIC d 

137+| 9 

15211 

1521® 



/abs,0.3-8keV ACIS d 




1081“ 

80l 2 ^ 

fPL d 

unabs, 0.5—8ke V 

12 ± 2 

9-9lo.2 

10 ± 1 

S-Slls 

H+5 

xl -4 

^O.CskeV (!° 3 ° er § s- 1 ) 

1.73 ± 0.25 

1 qo+0.16 

x.oo-o 03 

1.47 ± 0.16 

1 99+0.31 
i.zz-o 27 

1 63 +0 ' 76 
-*-•0*5—0.54 

Abol.CBB (10 31 ergs' 1 ) 

19.3±g 2 i 7 

8-4±l;5 

8 -611:2 

4 q+t-0 

0.7 

39.91 292 

Abol.HBB (10 3 ° erg s- 1 ) 

4.9±^J 

q o+i-0 
—0.6 

q 9 + 1-0 
O.Z_ 0 ? 

1 4 +0 - 7 

X - 4 —0.5 

3.6l®i 

X 2 /dof 

0.87 

0.99 

0.96 

0.96 

0.79 

dof 

87 

221 

274 

274 

52 

p n 


1 (fixed) 

1 (fixed) 



-Emosi 


1.06 ± 0.02 

1.06 ± 0.02 



AmOS2 


1.04 ± 0.02 

1.04 ± 0.02 



Aacis 




1 (fixed) 



a A model of two BBs and a PL was used. Uncertainties indicate the 90% confidence level. 
Luminosity uncertainties do not include the distance uncertainties. Luminosity and radius 
calculations were done using a distance value of 350 pc and the exact best-fit parameters before 
rounding. 

b The following best-fit parameters are given: absorbing hydrogen column density Ah, BB tem¬ 
perature kT , normalization A/pl of the PL component in units of 10“ 5 photons keV^ 1 cm -2 s —1 
at 1 keV, radius R for the equivalent sphere of the BB emission, photon index V of the PL com¬ 
ponent, fluxes /, bolometric luminosities of the BBs, the luminosity of the PL component in 
the energy range 0.5 — 8 keV, reduced x 2 , the number of degrees of freedom, and the calibration 
factors T in the case of tied fits, 

c pn, MOS1, MOS2, and ACIS were fitted together. The Ah of Chandra was tied to the XMM 
N#, and the XMM photon index was tied to the Chandra photon index. 
d Fluxes are given in units of 10 —14 erg cm -2 s —1 . 


Unfolded Spectrum 



Energy (keV) 

gas5167 30-Mar-2015 14:14 

Fig. 5.— The unfolded spectrum of the pulsar for the pn (black), 

MOS1 and MOS2 (combined, red) and ACIS (blue) instruments us¬ 
ing the spectral model fit where A# and the PL photon indices are 
tied for the XMM-Newton and Chandra data (columns 4 and 5 in 
Table [4|l. Note that we do a rigorous simultaneous fit of all individ¬ 
ual 7 spectra, but for better visibility, we show here the combined 
events of the pn, MOS and ACIS instruments, respectively. 



Fig. 6 .— Blackbody temperature confidence contours (68%, 
90%, 99 %) for the combined XMM-Newton/Chandra fit (A# and 
the PL photon indices are tied, see columns 4 and 5 in Table 0. 
The XMM-Newton contours are black, Chandra contours are red. 


emission. 


Since the work by iDe Luca et al.l (|2005l ) 






















new NS atmosphere models have become available, 
e.g. , the NS M a gnetic Atmosphe r e Mod e ls (NSMAX G) 
by IMori fc Hoi (12001 : IHo et, all (120081) : iHol (IM1 as 
implemented in XSPEC. We used the pn data to carry 
out fits using a NSMAXG+PL for a magnetic field 
B = 10 12 G, fixed NS mass of 1.4M 0 , NS radius of 
10 km and distance of 0.35 kpc. We checked fits with a 
hydrogen atmosphere (specfiles 1200, 123100, 123190 in 
XSPEC), carbon atmosphere (specfile 12006), oxygen 
atmosphere (specfile 12008), neon atmosphere (specfile 
12010). The fits were generally not acceptable (reduced 
X 2 > 2 for 89 dof) even if we allowed the mass, radius 
and distance to vary. We obtained the best fit for the 
hydrogen atmosphere (specfile 1200). However, the 
reduced y 2 = 1.4 for 89 dof is significantly worse than 
the one of the 2BB+PL fit (\ 2 = 0.9 for 87 dof), and 
there are substantial systematics in the soft part of the 
spectrum. Furthermore, the obtained normalization, 
J\T = f?Em/^NS = 34 (where i?Em indicates the size of 
the emission region), is rather unrealistic. Instead of 
a NSMAXG+PL we also tried NSMAXG+NSMAXG 
for this hydrogen atmosphere to consider possible hot 
spots, but the fit was not accetable (y 2 = 1.6 for 89 
dof, systematic residuals) either. We further checked 
a three-component model, NSMAXG+NSMAXG+PL 
for the hydrogen atmosphere (specfile 1200), but the 
fit was not accetable (y 2 = 1.3 for 87 dof, systematic 
residuals, Af = = 395l® 3s ) either. Overall, 

for a two- or three-component model approach, the 
NSMAXG (hydrogen) models seem not to be applicable 
to PSRB1055—52. 

4. DISCUSSION 

4.1. Extended Emission 

We detected extended emission around PSRB1055—52 
with a significance of 4 a (Table [2]) in the SEB for a region 
4'.'9 — 20". In Figure O the presence of the radial surface 
brightness enhancement of the observational data with 
respect to the MARX simulation is already visible at 
smaller radii, starting from T/5. Based on Figure [3l we 
did not use these smaller radii in order to optimize the 
significance analysis. However, the investigation of the 
V! 5 — 4" region around the pulsar (Figure [2] Table [3]) 
resulted in finding of a 3cr asymmetry in the present 
emission (there is no significant asymmetry at larger 
radii), with most counts in the SEB as well. 

It is difficult to derive further spectral and spatial 
information for the extended emission from this ACIS 
observation, because even in the BEB there are only 
« 30 (background-subtracted) counts for the chosen 
extraction region listed in Table [21 Since the extended 
X-ray emission at 4'.'9 — 20" is only significant in the 
SEB, a possible explanation for the excess could be a 
dust scattering halo, whose spectrum is generally softer 
than the central source. So far, no such halos have 
been reported for any of the neutron stars which have a 
similar low interstellar absorption (A# ~ 10 20 cm -2 ). 

As described in Section [3+1 we did not find any sig¬ 
nificant flux enhancement for angular separations > 20". 
Since detected dust scattering halo profiles are usually 


very broad (see, e.g.. iPredelil fe Schmittl [l~995h . it is in¬ 
teresting to check whether the surface brightness profile 
around PSRB1055—52 is consistent with adust scatter¬ 
ing halo profile. Based on th e works by iDraind ([2003T I 
and iPredehl fc Schmitt] (119951) , we use the following for¬ 
mula for the halo flux, i 7 haio(0i, # 2 , E), in an annulus 
9 1 < 6 h < 62 at a photon energy E: 

F halo (9 1 ,9 2 ,E) = E PSR (E)r sca (lkeV)^=^, (1) 

where Epsr is the flux of the pulsar, r SC a(l keV) is the 
scattering optical width at 1 keV, and the constant 0 de¬ 
pends on the dust model (we use 0 = 36 0 " as derived b y 
Drainel 120031 for a dust model bv IWeingartner fe Drain el 
200111 . Formula dT]) is only applicable for small scattering 
angles, 02 < 0- In order to use the measured count 
numbers instead of fluxes, one has to convolve the 
respective fluxes with the instrument response. The 
transformed formula for count numbers is very similar 
to equation ([!]) (for details, see, e.g., Birzan et al. 2015, 
submitted), but instead of (1/E) uses (1/E) e ff, the 
mean inverse energy of the central source photons. For 
PSRB1055—52, (1/E) eff = 1.3keV -1 . 

From the halo flux in one annulus, one can now 
estimate the expected count enhancement in another 
annulus. Because of the CCD gaps, we estimate the ex¬ 
pected count enhancement in a fan region, e.g., from 20" 
to 40" (see Fig. Q] where a fan region up to 60" is shown). 
From the enhancement of 31 counts in the BEB and 
27counts in the SEB in the respective annuli (Tabled, 
we estimate expected enhancements of 27 ± 5 counts 
(BEB) and 13 ± 4 counts (SEB) for the fan region. 
From our MARX simulation and the background rate 
measurement, we expect that 94 ± 11 counts (BEB) and 
10 ± 3 counts (SEB) actually come from the background 
and, negligibly, from the PSF wings. Hence, one would 
expect a total count number of 121 ± 12 counts (BEB) 
and 23 ± 5 counts (SEB). In the respective fan region of 
the actual data, however, we find 98 ± 10 counts (BEB) 
and 15 ± 4counts (SEB), i.e., the measured difference 
between halo prediction and data is 23 ± 16 counts 
(BEB) and 8 ± 6counts (SEB). Though we find notably 
fewer counts than one would expect in presence of a 
dust halo, the strong background contribution renders 
the significance of the difference between prediction and 
actual data to < 2 a in both energy bands. Hence, a 
dust halo still cannot be excluded. We note that such 
profile estimates are obviously biased with respect to 
the assumed dust grain properties and dust distribution 
along the line of sight. A sharper surface brightness 
profile could be produced, for example, by a localised 
dust clump very close to the pulsar. However, there 
is no obvious indication for such a dust clump in the 
direction of PSRB105 5—52 in 3D ISM models (see, e.g., 
lLallement et al.ll2014f) . 

As shown by IPredehl k. Schmitt ([1995( 1. one 
can assume r sca (l keV) = SAh, 22 with Ah, 22 = 

Ah/(10 22 cm -2 ) and S a constant on t he order of 1 
(e.g., S ~ 0.5; IPredehl fe Schmittl 119951 ). Aiming to 
constrain S, we use the measured BEB pulsar counts, 
1858 + 43 (Section[3|), the determined BEB counts of the 
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enhancement in the annulus 4.88" — 13.10", 31.2 ± 9.4 
from Table [H and IVh, 22 = 0.034 ± 0.001 from the 
Chandra/XMM fit of the pulsar spectrum. We obtain 
S = 10.5 ± 3.2. This value appears to be substantially 
larger than one expects for a dust scattering halo. 
Due to the low count number in the enhancement and 
the resulting large 3cr error of S , however, we cannot 
entirely exclude that a dust scattering halo is the only 
or main source of the found flux enhancement around 
PSRB1055—52. 

The 3cr asymmetry of X-ray count distribution in 
the 1"5 — 4" region around the pulsar could be ei¬ 
ther a statistical noise feature or indication of a very 
faint and very compact PWN. Since the slight count 
enhancement is located in the direction of the proper 
motion, one could su spect a b ow-s hock PWN. From 
optical observations, IMignani et all (I2010t ) measured 
the proper motion of PSRB1055—52 as p = 42 ± 5 
mas yr -1 . This translates into a transverse velocity 
v± = (70 ± 8)kms -1 at a distance of 350 pc. As this 
velocity is considerably larger than the typical sound 
speed in the ISM, cs ~ 10 — 30kms -1 , one would 
indeed expect that the pulsar should create a bow-shock 
PWN. Using the transverse velocity as a lower limit 
to the total velocity, a typical ambient density of 
the ISM, pism ~ 2 x 10 -24 gcm -3 , and ap plying the 
stand-off radius R s formula (e.g., IChatteriee fe Cordesl 
I2002T) . R 2 S = -E (47rcpisM'^totai) -1 ’ one can estimate 
R s < 2.7 x 10 16 cm, corresponding to < 5" at d = 350 pc 
for the expected stand-off radius of a bow-shock PWN 
around PSRB1055—52. Here, smaller values of R s 
would be expected if Utotai > u_l, he., if a significant 
radial velocity component is present. Overall, the found 
asymmetry - if real - is consistent with, or slightly 
smaller than the expected bow shock dimensions. 
However, given the faintness of this asymmetry emission 
and the preponderance of soft counts, background 
fluctuations are equally likely. 

One could speculate about a model which could 
accommodate the different radial extents by assuming 
that PSRB1055—52 has a significant radial velocity, 
i.e., the pulsar is moving towards us or away from us. 
From simulations of bow-shock nebulae, it was found 
that the radial extent of the cylindrical PWN t ail can 
be es timated as Rt a u ~ 4 R s (e.g., IBucciantini et al.l 
l2005f) . From these models, one expects the shocked 
pulsar wind to flow with nearly relativistic velocities. 
Therefore, considering beaming and projection effects 
for a pulsar moving along the line of sight, one would 
expect to see a PWN with an extent of ~ 4 R s if the 
pulsar moves away from us and an extent of ~ R s if the 
pulsar moves toward us. If the pulsar moves away from 
us with Utotai u_l, we would expect a nearly round 
PWN. In the case of PSRB1055—52, additional radial 
movement away from us could explain the observed 
features: the 4"9 — 13" flux enhancement and the small 
1"5 — 4" asymmetry. 

The estimated unabsorbed BEB flux of the enhance¬ 
ment (see Section m , corresponds to a luminosity of 
about TggB ~ 1 x 10 29 d 3 50 ergs -1 or ~ 2 x 


10 29 d 2 50 ergs -1 with ^350 = d/350 pc for assumed pho¬ 
ton indices of 2 or 4, respectively. If the extended emis¬ 
sion is mainly due to a PWN, the corresponding effi¬ 
ciency ))pwn = T pwn /E is between 3 — 8 x 10 -6 . This 
efficiency is lower than for typical PWN around young 
pulsars, but in fact very simil ar to t he value of Geminga’s 
PWN (e.g., Figure 5 bv lKargaltsev fc Pavlovl[2n08l ). The 
found prevalence of soft counts could be either due to an 
unusually soft PWN spectrum or just due to the faint¬ 
ness of the PWN (for a typical power law slope, more 
flux contribution is expected from soft than from hard 
photons). The current count numbers are too low to rea¬ 
sonably constrain the spectrum of the flux enhancement 
and disentangle putative PWN and scattering dust halo 
contributions. As explained above, a major contribution 
from a dust scattering halo seems unlikely, but cannot 
entirely be excluded. 

4.2. Pulsar 

An obvious reason for the differences between the 
XMM-Newton EPIC and Chandra ACIS results for 
the pulsar spectrum could be cross-calibration effects. 
Cross-calibration between A CIS a n d EP IC has been in- 
yestig a ted bv ISch ellenber ger et al.l (120151) : IK ettula e t al.1 
(|2013lklTsunmoto et al.1 (l2011l) : INevalainen et all (120101. 
These studies compare, e.g., spectral fit results for ab¬ 
sorbed (Njj ~ a few 10 2 O cm -2 ) relatively hard sources, 
usually galaxy clusters. They consistently find higher 
temperatures for ACIS than for EPIC, in particular 
for EPIC-pn. In th eir soft X-ray band (0.2 — 2.0 keV), 
IKettula et all (I2013H reported 18% higher ACIS (APEC 
model) temperatures than obtained with EPIC-p n. We 
refer to Figures 5, 6 and 11 by ISchellenberger et aH 
(12015T) for an illustration of the higher- ACIS-temperat ure 
effect. In addition, Figure 13 by ISchellenberger et~akl 
(120151 ) shows that ACIS results usually imply a larger 
Nh than pn results. For soft X-ray sources, there are 
no comparable systematic studies to our knowledge. An 
example for a comparison between soft spectra from 
ACIS and EPIC (with sufficient number of counts) is 
t he stu dy of the isolated neutron star RBS 1223 by 
lHaberl et al.l (120031 ) . In that work, the ACIS-80 data 
from 2000 and the EPIC data from 2001 and 2003 
were fit with an absorbed blackbody with a broad 
gaussian absorption line. The BB temperature was only 
slightly higher for ACIS-S (87.8 ± 1.0 versus 85.8 ± 0.5 
for EPIC), and the inferred Njj limit was smaller 
(< 1.6 x 10 2 ° cm -2 versus < 4.1 ± 0.1 x 10 2 ° cm -2 
for EPIC-pn) than the EPIC values. Hence, for these 
particular observations of a mildly absorbed soft source, 
one obtains consistent BB temperatures, and consistent 
or smaller Nh values for ACIS with respect to EPIC 

in contrast to the findings in the systematic studies 
of the galaxy clusters. Without a systematic study of 
mildly absorbed soft sources it is therefore difficult to 
accurately assess the actual cross-calibration effects for 
a comparison of ACIS / EPIC spectral fit results. 

For RBS1223, lHaberl et al.1 (120031 ) found a difference 
of about 10 % in t he absorbed flux (i n the same direction 
as our results). iNevalainen et al.1 (120101) (their Figure 

8 Note that there are slig ht differences for ACIS-S and ACIS-I 
comparisons (e.g., Fig. 2 by Schellcnbergcr ct al. 2015) 
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20, Table 11) reported that flux differences between 
ACIS and pn results were only 2% in their soft band 
(0.2 — 2.0 keV) using ACIS data obtained before 2008. 
For PSRB1055—52, the flux difference between the 
EPIC and ACIS results is 29% and 47% for tied (Nh, 
PL photon index) and free parameter fits, respectively 
(Table |U). Even considering errors, this appears to be 
much larger than a reasonable cross-calibration effect 
expected from previous analyses. Our Chandra obser¬ 
vations helped to evaluate the effect of contaminating 
neighbor sources on the flux. Inclusion of Source A 
in the ACIS-only fit of the 2BB+PL model resulted 
in an observed flux of 8.ljj'} x 10 -13 ergcm -2 s -1 in 
comparison to 8 . 0 l 2 '! x 10 -13 ergcm -2 s -1 if this source 
was not included. Thus, one would expect the XMM- 
Newton flux to be higher by only 1% (with a formal 
error of 37%) if Source A was not subtracted (it was for 
MOS1/2). The formal errors are large due to uncertain 
spectral fit constraints. However, the corresponding 
increase in the count numbers of 4 ± 3% (Source A has 
77 counts, the pulsar has 1912; see Section EO) support 
the finding that less than 10 % of the flux difference 
can be attributed to inclusion of Source A. The other 
possible contaminating sources would have an even 
smaller (factor of 4) effect on the XMM-Newton results. 
Hence, contaminating neighbor sources and a formal 
cross-calibration error of 10 % cannot account for an 
overall flux difference of 29% or more. 

With respect to the cross-calibration studies men¬ 
tioned above, one needs to consider possible calibration 
changes over time, too. Whil e EPIC-pn w as found to be 
very stable over time (iSartore et al.ll2012l) . this is not the 
case for ACIS. A contaminant has been accumulating 
on the optical-blocking filters of the ACIS detectors 
with a different rate over time. While the effect on the 
inferred spectral parameters of an X-ray source should 
be minimized by the contamination model implemented 
in the ACIS data reduction, an unavoidable uncertainty 
in the contamination correction can lead to a different 
measured flux and can offset the derived spectral param¬ 
eters. A underestimated contamination, for example, 
would lead to a larger-tlian-real Nh in the spectral fit 
results. Recently, we discussed a 4%-6% decrease of 
the absorbed flux (and apparent slight Nh increase) for 
the Central Compact Object in the Cas A supernova 
remnant whic h was measured with ACIS-I between 2006 
and 2012 ( Posselt et al. 2013 ). We showed that, in 
principle, the uncertainty in the ACIS calibration in a 
small energy range could be responsible for this small 
effect. In the case of PSRB1055—52, it is intriguing 
that the temperatures and emission area sizes of the 
ACIS-only fit (free and very large Nh which could be 
partly due to the ACIS contamination) are much closer 
to the EPIC values than the numbers obtained if Nh 
is tied to the much smaller EPIC value. Currently, 
we cannot explain the surprisingly large flux difference 
between the ACIS and EPIC results. We suspect a 
cross-calibration issue, but this can only be probed 
with further X-ray observations, preferably with the 
time-stable EPIC-pn. 

Our Chandra data are also useful to assess the effect 


of contaminating neighbor sources on the measured pho¬ 
ton index. As described in Section 13.21 exclusion of the 
neighbor source A results in a slightly softer PL. In par¬ 
ticular, using more restricted source extraction regions 
for the XMM-Newton data together with the ACIS data, 
we fin d a slightly larger photon index than lDe Luca et akl 
(120051) (1.9 ± 0.2 versus 1.7 ± 0.1) The PL component 
of the optical spectrum of PSR B1055—52 has a photon 
index of 2.05 ± 0.34 (iMignani et al.ll2010H . Though for¬ 
mally both components have similar photon indices, the 
new X-ray photon index does not change the fact that 
the extension of the X-ray PL component still overshoots 
the optical fluxes by a factor ~ 4. This is, however, sim- 
ilarjo most o ther pulsars detected in the optical (e.g., 
iMignanil 120111) . We note that at Fermi 7 -r ay ener g ies, 
the p hoton index is harder, T 7 = 1.1 ± 0.1 (|Abdo et alJ 
120131 : see also Table [T]) , indicating an additional spectral 
break between the X-ray and the y-ray spectral ranges. 

4.3. Comparison with the PWN properties of other 
pulsars 

Of the 13 nearby (< 2kpc) middle-aged rotation- 
powered pulsars that have sensitive public X-ray 
observations with Chandra or XMM-Newton (Table [lj , 
only 6 have a prominent (relatively bright, extended) 
X-ray PWN, while 7 others have none or very compact 
faint ones. The question arises what is the defining 
physical factor for the presence of a prominent X-ray 
PWN ? One needs to consider the intrinsic pulsar 
properties, exterior environment factors, as well as 
sight-line dependent geometric effects. Our selection for 
a restricted E range ensures that the spin-down energy 
alone cannot be the defining factor. PSR J1741—2054, 
for example, has a factor of 3 smaller E than the Three 
Musketeers, but a prominent X-ray PWN. Since the 13 
pulsars have very similar spin-down properties, their 
inferred magnetic dipole fields have similar strengths. 
Thus, one would expect similar magnetospheres, capable 
of producing similar amount of accelerated charged 
particles. The cutoff energy in y-rays, another indicator 
of the availability of sufficiently accelerated particles in 
the magnetosphere, does not show a correlation with 
(non-)detections of PWNe either (Table [0. 

The efficiency of particle acceleration in the mag¬ 
netosphere can depend on the angle between the 
rotational and magnetic axis, the obliquity a (e.g., 
iPliilionov et al.ll2014t iTchekhovskov et al.N2013l) . While 
axis align ment in pulsars has been suggested on Myr 
timescales (lYoung et al.ll 20 ltil) . the pulsars in Table Q] are 
too young for this effect to be relevant. Different axis 
orientations can produce different pattern of detections 
in radio, X-rays and y-rays (see, e.g., lAbdo et~aT1120131 
for a discussion about different radio-loud/radio-quiet 
fractions of Fermi pulsars with different E). For the 
considered middle-aged pulsars, there is no obvious 
correlation between the PWN (non-)detections and 
loudness or quietness in the radio and y-ray frequency 
ranges. Given the dependence of radio and y-ray 
detections on the respective beam widths and the 
angles between the direction of sight and the rotation 
axis, £, th e disentanglement of g eometric factors is 
non-trivial. iPierbattista et ahl (|2015l f used four different 
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magnetosphere models (PC - Polar Cap, SG - Slot 
Gap, OG - Outer Gap, OPC - One Pole Caustic) to fit 
radio and 7 -ray pulse profiles separately and together. 
Generally, all the models perform poorly, in particular if 
both wavelength ranges are considered together. With 
respect to PSRB1055—52 (very faint/compact PWN), 
we restrict the discussion on the influence of geometric 
factors to the exemplary comparison with the Geminga 
pulsar (very prominent PWN). 

For the Geminga p u lsar, t he best fitting model used 
by iPierbattista et al] ( 2015 1. the SG model, resulted 
in a = 42°, £ = 51°. Even this model, however, 
poorly reproduces G e minga ’s 7 -ray interpulse. Previ¬ 
ously, IWatters et al.1 (|2009f) also obtained constraints 
on the geometry of the Geminga pulsar, but did not 
show the respective fitted pulse profiles. Using the 
OG-model they derived a = 10° — 25°, £ = 85°; 
using a two-pole caustic magnetosphere model they 
obtained two solutions: a = 30 ° — 8 0°, £ = 90° and 
a = 90°, £ = 55° — 80°. iMalovl (Il998i l argued that the 
radio-quietness of Geminga indicates an aligned rotator. 
While the current constraints on Geminga’s a from its 
7 -ray pulse profile remain uncertain, the models seem 
to agree that £ > 50°. If one assumes that the direction 
of proper motion and rotation axis are crudely aligned, 
this constraint on £ means that Geminga’s total velocity 
is close to its measured transverse velocity. Geminga’s 
PW N consists o f two large lateral tails and an axial tail 
(e.g.- lPavlov et al.ll2010f ). The emission properties of the 
lateral tails resemble more those of jets than of a bow 
shock front, a fact which - together with the axial tail 
- seems to be indicative of a large a (Posselt et al., in 
preparation). 

For PSR B1 055—52 iLvne &; Manchester! (I1988T ) and 
iWeltevrede fc Wrightl (|2009t ) derive very similar obliq¬ 
uity values, a = 75°, £ = 67° — 69°, from radio 
polarization measurements. Assuming that rotation 
axis and sight line are aligned, iMalov fe Nikitinal (j201lH 
also used radio polarization measurements to obtain 
a = 8 ° for PSRB1055—52. Among the models used by 
IPierbattista et all (120151! to describe the radio and 7 -ray 
light curves of PSRB1055—52, the OG-model and the 
PC-model gave the best fit results with a = 77° , £ = 87° 
and a = 10°, £ = 7°, respectively. The radio pulse 
shape is poorly reproduced in the fits for both models. 
Overall, assuming an alignment of rotation and proper 
motion axes, the current constraints on PSRB1055—52 
seem to support either an oblique rotator, moving nearly 
perpendicular to the line of sight, or a nearly aligned 
rotator which moves along the line of sight. 

The size of a bow shock created by the PWN around 
a pulsar depends not only on the pulsar energetics, but 
also on the pulsar velocity and the density of the ISM 
(as we discussed for PSR B1055—52 in Section B~T1) . The 
transverse velocity of PSRB1055—52, 70kms _1 , is low in 
comparison to the mean value of other normal pulsars, 
246 ± 22kms~ 1 (IHobbs et al.1120031 . PSR J0659+1414 
which also does not have a prominent PWN, has a similar 
low transverse velocity, while Geminga’s t ransverse ve¬ 
locity is close to the mean value reported bv IHobbs et al.1 
(f2005il (see Table \Q. However, PSR J0358+5413 has a 


transverse velocit y of only SSkms -1 , y et also powers a 
prominent PWN (IMcGowan et al.ll2006T) . Clearly, differ¬ 
ent velocities alone cannot explain the diverse picture of 
PWN detections around middle-aged pulsars. 

The influence of the ISM density can be crudely 
judged by comparing Njj values. For instance, different 
Nh for similarly distant pulsars in the same sky area 
can indicate different ISM densities close to the pulsars 
themselves. The Nh values of the pulsars in Table Q] are 
all very similar, but of course local inhomogeneities in 
the ISM density may exist. It is particularly interesting 
to compare Geminga and PSR J0659+1414. These 
pulsars are very close to each other, and the local 3D 
map of the ISM does not reveal any substantial density 
variations between these two pulsar positions (at a 
spatial resolution of ~ 10 pc; lLallement et al.112014H . 
Though PSRB1055—52 is in a different part of the sky, 
its Nh value and the density in its local 3D map of the 
ISM resemble those of Geminga and PSR J0659+1414. 
The ISM conditions of the Three Musketeers appear to 
be similar enough that they cannot explain the pulsar’s 
different PWN manifestations. 

Overall, the governing factor for the different PWN 
properties of the middle-aged pulsars appears to be 
the obliquity of the pulsar in combination with the 
inclination to the observer. With respect to the 
Three Musketeers, we suggest that PSRB1055—52 and 
PSR J0659+1414 are nearly aligned rotators whose main 
velocity component is along the line of sight. In con¬ 
trast, Geminga appears to be an highly oblique rotator 
which moves nearly perpendicular to the line of sight. 
Further improve ments in pulsar magnetosp here models, 
e.g., those used bv IPierbattista et al.1 (|2015l) . and the cor¬ 
responding detailed analyses of 7 -ray, X-ray and radio 
pulse shapes are needed to confirm this hypothesis. 

5. CONCLUSIONS 

The vicinity of PSRB1055—52 does not show a 
prominent X-ray PWN. A 4<r flux enhancement has 
been detected in the SEB in an annulus 4'/9 — 20" 
around the pulsar. Taking also into account the X- 
ray count distribution at larger separations from the 
pulsar and properties of known dust scattering halos, 
a PWN interpretation seems more grounded than a 
dust scattering halo, though we cannot entirely exclude 
a contribution from the latter due to the low count 
numbers. We detected a slight (3cr) asymmetry in the 
count distribution l'/5 — 4" around the pulsar. This 
asymmetry is most pronounced in the BEB and located 
in the direction of the pulsar proper motion. While this 
emission is consistent with a bow shock PWN compo¬ 
nent due to the pulsar’s transverse velocity, we cannot 
rule out background count fluctuations or MARX model 
uncertainties as other causes for the detected asymmetry. 

Comparing the properties of nearby middle-aged 
pulsars, we identified differences in geometrical axis 
orientations as the most likely reason for the fact that 
some of these pulsars have very prominent X-ray PWNe 
while their - otherwise very similar - cousins have 
none or very faint/compact PWNe. The faint extended 
emission around PSRB1055—52 is consistent with a bow 
shock around an aligned rotator moving nearly along 
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the line of sight. 

For PSRB1055—52, we found a 30% (or larger) flux 
decrease between the 2000 XMM-Newton and the 2012 
Chandra observation which could formally be described 
by smaller (but hotter) thermal emission areas. We 
suspect a cross-calibration issue, but this needs to be 
proven with further X-ray observations, preferably with 
the time-stable XMM -Newton EPIC-pn. 


The scientific results reported in this article are based 
on observations made by the Chandra X-ray Observa¬ 
tory. Support for this work was provided by the National 
Aeronautics and Space Administration through Chan¬ 
dra Award Number G02-13090X issued by the Chan¬ 
dra X-ray Observatory Center, which is operated by the 
Smithsonian Astrophysical Observatory for and on be¬ 
half of the National Aeronautics Space Administration 
under contract NAS8-03060. 
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